A metagenomic investigation of phytoplasma diversity in Australian vegetable growing regions.

In this study, metagenomic sequence data was used to investigate the phytoplasma taxonomic diversity in vegetable-growing regions across Australia. Metagenomic sequencing was performed on 195 phytoplasma-positive samples, originating either from historic collections (n=46) or during collection efforts between January 2015 and June 2022 (n=149). The sampled hosts were classified as crop (n=155), weed (n=24), ornamental (n=7), native plant (n=6), and insect (n=3) species. Most samples came from Queensland (n=78), followed by Western Australia (n=46), the Northern Territory (n=32), New South Wales (n=17), and Victoria (n=10). Of the 195 draft phytoplasma genomes, 178 met our genome criteria for comparison using an average nucleotide identity approach. Ten distinct phytoplasma species were identified and could be classified within the 16SrII, 16SrXII (PCR only), 16SrXXV, and 16SrXXXVIII phytoplasma groups, which have all previously been recorded in Australia. The most commonly detected phytoplasma taxa in this study were species and subspecies classified within the 16SrII group (n=153), followed by strains within the 16SrXXXVIII group ('Ca. Phytoplasma stylosanthis'; n=6). Several geographic- and host-range expansions were reported, as well as mixed phytoplasma infections of 16SrII taxa and 'Ca. Phytoplasma stylosanthis'. Additionally, six previously unrecorded 16SrII taxa were identified, including five putative subspecies of 'Ca. Phytoplasma australasiaticum' and a new putative 16SrII species. PCR and sequencing of the 16S rRNA gene was a suitable triage tool for preliminary phytoplasma detection. Metagenomic sequencing, however, allowed for higher-resolution identification of the phytoplasmas, including mixed infections, than was afforded by only direct Sanger sequencing of the 16S rRNA gene. Since the metagenomic approach theoretically obtains sequences of all organisms in a sample, this approach was useful to confirm the host family, genus, and/or species. In addition to improving our understanding of the phytoplasma species that affect crop production in Australia, the study also significantly expands the genomic sequence data available in public sequence repositories to contribute to phytoplasma molecular epidemiology studies, revision of taxonomy, and improved diagnostics.

The observation of taxonomic boundaries for the 16SrII and 16SrXXV phytoplasmas using genome-based delimitation.

Within the 16SrII phytoplasma group, subgroups A-X have been classified based on restriction fragment length polymorphism of their 16S rRNA gene, and two species have been described, namely 'Candidatus Phytoplasma aurantifolia' and 'Ca. Phytoplasma australasia'. Strains of 16SrII phytoplasmas are detected across a broad geographic range within Africa, Asia, Australia, Europe and North and South America. Historically, all members of the 16SrII group share ≥97.5 % nucleotide sequence identity of their 16S rRNA gene. In this study, we used whole genome sequences to identify the species boundaries within the 16SrII group. Whole genome analyses were done using 42 phytoplasma strains classified into seven 16SrII subgroups, five 16SrII taxa without official 16Sr subgroup classifications, and one 16SrXXV-A phytoplasma strain used as an outgroup taxon. Based on phylogenomic analyses as well as whole genome average nucleotide and average amino acid identity (ANI and AAI), eight distinct 16SrII taxa equivalent to species were identified, six of which are novel descriptions. Strains within the same species had ANI and AAI values of >97 %, and shared ≥80 % of their genomic segments based on the ANI analysis. Species also had distinct biological and/or ecological features. A 16SrII subgroup often represented a distinct species, e.g., the 16SrII-B subgroup members. Members classified within the 16SrII-A, 16SrII-D, and 16SrII-V subgroups as well as strains classified as sweet potato little leaf phytoplasmas fulfilled criteria to be included as members of a single species, but with subspecies-level relationships with each other. The 16SrXXV-A taxon was also described as a novel phytoplasma species and, based on criteria used for other bacterial families, provided evidence that it could be classified as a distinct genus from the 16SrII phytoplasmas. As more phytoplasma genome sequences become available, the classification system of these bacteria can be further refined at the genus, species, and subspecies taxonomic ranks.

Iodixanol density gradients as an effective phytoplasma enrichment approach to improve genome sequencing.

Obtaining complete phytoplasma genomes is difficult due to the lack of a culture system for these bacteria. To improve genome assembly, a non-ionic, low- and iso-osmotic iodixanol (Optiprep™) density gradient centrifugation method was developed to enrich for phytoplasma cells and deplete plant host tissues prior to deoxyribonucleic acid (DNA) extraction and high-throughput sequencing (HTS). After density gradient enrichment, potato infected with a 'Candidatus Phytoplasma australasia'-related strain showed a ∼14-fold increase in phytoplasma HTS reads, with a ∼1.7-fold decrease in host genomic reads compared to the DNA extracted from the same sample without density gradient centrifugation enrichment. Additionally, phytoplasma genome assemblies from libraries equalized to 5 million reads were, on average, ∼15,000 bp larger and more contiguous (N50 ∼14,800 bp larger) than assemblies from the DNA extracted from the infected potato without enrichment. The method was repeated on capsicum infected with Sweet Potato Little Leaf phytoplasma ('Ca. Phytoplasma australasia'-related strain) with a lower phytoplasma titer than the potato. In capsicum, ∼threefold more phytoplasma reads and ∼twofold less host genomic reads were obtained, with the genome assembly size and N50 values from libraries equalized to 3.4 million reads ∼137,000 and ∼4,000 bp larger, respectively, compared to the DNA extracted from infected capsicum without enrichment. Phytoplasmas from potato and capsicum were both enriched at a density of 1.049-1.058 g/ml. Finally, we present two highly contiguous 'Ca. Phytoplasma australasia' phytoplasma reference genomes sequenced from naturally infected Solanaceae hosts in Australia. Obtaining high-quality phytoplasma genomes from naturally infected hosts will improve insights into phytoplasma taxonomy, which will improve their detection and disease management.

'Candidatus Phytoplasma stylosanthis', a novel taxon with a diverse host range in Australia, characterised using multilocus sequence analysis of 16S rRNA, secA, tuf, and rp genes.

In Australia, Stylosanthes little leaf (StLL) phytoplasma has been detected in Stylosanthes scabra Vogel, Arachis pintoi Krapov, Saccharum officinarum L., Carica papaya L., Medicago sativa L., and Solanum tuberosum L. The 16S rRNA gene sequence of StLL phytoplasma strains from S. scabra, C. papaya, S. officinarum and S. tuberosum were compared and share 99.93-100 % nucleotide sequence identity. Phylogenetic comparisons between the 16S rRNA genes of StLL phytoplasma and other 'Candidatus Phytoplasma' species indicate that StLL represents a distinct phytoplasma lineage. It shares its most recent known ancestry with 'Ca. Phytoplasma luffae' (16SrVIII-A), with which it has 97.17-97.25 % nucleotide identity. In silico RFLP analysis of the 16S rRNA amplicon using iPhyClassifier indicate that StLL phytoplasmas have a unique pattern (similarity coefficient below 0.85) that is most similar to that of 'Ca. Phytoplasma luffae'. The unique in silico RFLP patterns were confirmed in vitro. Nucleotide sequences of genes that are more variable than the 16S rRNA gene, namely tuf (tu-elongation factor), secA (partial translocation gene), and the partial ribosomal protein (rp) gene operon (rps19-rpl22-rps3), produced phylogenetic trees with similar branching patterns to the 16S rRNA gene tree. Sequence comparisons between the StLL 16S rRNA spacer region confirmed previous reports of rrn interoperon sequence heterogeneity for StLL, where the spacer region of rrnB encodes a complete tRNA-Isoleucine gene and the rrnA spacer region does not. Together these results suggest that the Australian phytoplasma, StLL, is unique according to the International Organization for Mycoplasmology (IRPCM) recommendations. The novel taxon 'Ca. Phytoplasma stylosanthis' is proposed, with the most recent strain from a potato crop in Victoria, Australia, serving as the reference strain (deposited in the Victorian Plant Pathology Herbarium as VPRI 43683).